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Abstract 

We apply the averaging method to analyze spatio-temportal structures in nonlinear Schrodinger 
equations and thereby study the dynamics of quasi-one- dimensional collisionally inhomoge- 
neous Bose-Einstein condensates with the scattering length varying periodically in spatial 
and crossing zero. Infinitely many (positive measure set) modulated amplitude waves (pe- 
riodic and quasi-periodic), which are instable, can be proved to exist by adjusting the inte- 
gration constant c on some open interval. Finally, some numerical simulations support our 
results. 

Keywords: Modulated amplitude waves; Gross-Pitaevskii equations; Collisionally 
inhomogeneous Bose-Einstein condensates; Averaging method 
PACS: : 05.45.-a, 03.75.Lm, 05.30.Jp, 05.45.Ac 



1. Introduction 

Since the experimental realization of Bose-Einstein condensates (BECs) in the mid-1990s 
[H0], the study of matter- wave patterns includingexistence and stability in BECs has drawn 
a great deal of interest from experimentalists and theorists . 

In atomic physics, the Feshbach resonance of the scattering length of interatomic inter- 
actions is used for control of Bose-Einstein condensates j9[ [Io| . We consider the main model 
of this paper for Feshbach resonance given by the perturbed Gross-Pitaevskii (GP) equation 



of the dimensionless form 11 



= + s( x M 2 ^ + H*)*, (1.1) 

where the nonlinearity coefficient g(x) varies in space. In Eq. (11. ip . ip is the mean- field 
condensate wave function (with density \ip\ 2 measured in units of the peak ID density n ), 
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x and t are normalized, respectively, to the healing length £ = H/ ^/no\gi\m and £/c (where 
c = ^--\/^o|fi , i|/ m is the Bogoliubov speed of sound), and energy is measured in units of the 
chemical potential S = gin . In the above expressions, g\ = 2frjjj±a , where u± denotes the 
confining frequency in the transverse direction, and ao is a characteristic (constant) value of 
the scattering length relatively close to the Feshbach resonance. Finally, V(x) is the rescaled 
external trapping potential, and the x-dependent nonlinearity is given by g(x) = a(x)/a , 
where a(x) is the spatially varying scattering length. 

In the past few years, GP equation (II. IB h as been widely studied, such as the stability and 
dynamics of bright, dark solitary waves |12l-ll5j and modulated amplitude waves (MAWs) 
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In order to study the dynamics of BECs with scattering length subjected to a spatially 



periodic variation, Porter and Kevrekidis et al. 11] transform equation (II. ip into a new GP 



equation with a constant coefficient and an additional effective potential 



V eff{X 



dt 2 dx 2 

1 f" (/') 2 , /' d 



2 / p f dx 



with f(x) = ^/g(x), then the transformed equation was investigated. For weak underlying 
inhomogeneity, the effective potential takes a form resembling a superlattice, and the am- 
plitude dynamics of the solutions of the constant-coefficient GP equation obey a nonlinear 
generalization of the Ince equation. In the small-amplitude limit, they use averaging to 
construct analytical solutions for modulated amplitude waves (MAWs), whose stability was 
subsequently examined using both numerical simulations of the original GP equation and 
fixed-point computations with the MAWs as numerically exact solutions. However, men- 
tioned in their paper, the transformation z/5 = ^g(x)tp applies only in the case when g(x) 
does not cross zero. A natural question is that for general periodic function g(x), whether 
the similar results upon the dynamics can be obtained. On the other hand, the phases 



of MAWs considered in [11[ are trivial, which are corresponding to standing waves. Thus, 
another question is that whether the MAWs with nontrivial phases can exist. 

With these questions discussed above, in this paper we investigate the existence and 
stability of MAWs with nontrivial phases in collisionally inhomogeneous BECs modeled 
by GP equation (II. ip for general small periodic function g(x). The method is based on 
averaging, and we use the averaging principle to replace a GP equation by the corresponding 
averaged system. Along this paper, we assume that g(x) and V(x) are analytic and periodic 
functions with the least positive period 

The rest paper is organized as follows. In Section [21 we introduce modulated amplitude 
waves involving periodic and quasi-periodic, and an averaging theorem is obtained in Section 
|3j In Section H] we investigate the existence and stability of equilibrium points for the aver- 
aged system and thereby study the periodic orbits and a numerical simulation is presented 
as prescribed parameters in Section |5j Finally, we summarize our results in Section [6j 
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2. Coherent structure 



we consider uniformly propagating coherent structures with the ansatz 

x/j(t, x) = R(x) exp(z[6(x) - fit}), (2.1) 

where R(x) e M gives the amplitude dynamics of the condensate wave function, 6(x) deter- 
mines the phase dynamics, and the "chemical potential" /i, defined as the energy it takes 
to add one more particle to the system, is proportional to the number of atoms trapped in 
the condensate. When the (temporally periodic) coherent structure (12.11) is also spatially 
periodic, it is called a modulated amplitude wave (MAW) (lil . Il7j |. Similarly, a solution 
of the equation (11.11) with the (temporally periodic) coherent structure (12. ip is called a 
quasi-periodic modulated amplitude wave (QMAW) if it is also spatially quasi-periodic. 

Inserting (12.11) into (II. ip . we obtain the following two couple nonlinear ordinary differ- 
ential equations 

r 2 

R" + 5R + eg(x)R 3 + eV{x)R = 0, (2.2) 

R d 

B" + 2Q'R'/R = => Q'(x) = (2.3) 

R z 

where 

eg(x) := —g(x), eV(x) := — V(x) 
and the integration constant c, determined by the velocity and number density, plays the 



role of "angular momentum" [18 



In case of c = 0, the phase of the condensate wave function (standing wave) is trivial and 
constant. In the general case, c ^ 0, the system (12.21) becomes more complicated and the 
phase is no longer constant [19|. Even the amplitude -R(x), a solution of (I2.2p . is T-periodic, 
the corresponding condensate wave function ip(x, t) may be not periodic, but quasi-periodic, 



with respect to the spatial variable x [20 



3. Averaging theorem 

Rewrite equation (I2.2p in the planar equivalent form 

R' = S 

c 2 (3 1) 

S' = -5R + — - eg(x)R 3 - eV{x)R. v ' ; 

Generally, averaging method involves two steps: transforming to standard form; solving the 
averaging equation. In order to proceed we need to transform (13. ip to a standard form for 
the method of averaging. 
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Lemma 3.1. Under the transformation f :Tx ( y T' +°°) ~~ * (0' +°°) x ^ defined by 



COS VOX + U) Sill VOX 



' cos" vox 



sm 2 (VSx + 

op 4 



system (13. ip changes into a new system 



p'=s 



9(x) 3 



c 2 c 2 
-(1 + — ) sin2(v / 5x + 0) + -(1 - — ) sin4(v / 5x 

V(x) sin2(v / 5x + ^ 



2^ 



0' = £ 



- (3 + cos 4( V5x + 9)) + 262p2{5p4 _ g2) cos 2( V5: 



(3.2) 



85 2 p 2 



z + 0) 



+ ^w( 1 + ^ cos2 (^ +( '))} 



i/ie new coordinates (9,p) in the half-plane T x (y y,+oo). 

The transformation \1/ arises from the variation of constant by using the solutions of the 
unperturbed system (e = 0), and p plays the role of "energy". When taking the integration 
constant c = 0, the transformation \1/ is the usual change of polar coordinates in the half 
plane. The proof the Lemma [3.11 follows from the basic computation (maybe lengthy), and 
it can be found in 20]. 

Now write the T-periodic functions g(x), V(x) as the Fourier series 



g{x) = go + ^^(ttfc sm2kV~5x + (3^ cos 2kVdx), (3.3) 

k=l 

oo 

V(x) = Vp + y^(afc sin2fcv / fa + 6 fc cos2kVSx). (3.4) 

k=l 

After inserting (13. 3p and (13 .4p into (13. 2 p and then multiplying the right side of (I3.2p by 1/T 
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and integrating from to T, we obtain the averaged system 

p ' = £ { %7lp~ A Sin(2 ^ + 0l) + 471 5 sin(2 ^ + 02) } := eFl ^ p) 

ff J Zgoi^ + c 2 ) (5 2 p 8 + c 4 ) (3.5) 
where 



0i = arctan (0 X = ^ • sign(0i), if ft = 0), 
Pi ^ 

0i = arctan — ^, (0i = — ■ sign(ai), if &i = 0). 

Theorem 3.1. [ Averaging theorem ] There exists a c r ,r > 2, change of variables 

p = p + ewi(8, p, x,e), 8 = 8 + ew 2 (8, p, x, e) 

with wi,W2 T -periodic functions of x, transforming (13.21) into 

p' = eFi (0, p) + e 2 gi (8, p, x,e) 
8> = eF 2 (8,p)+e 2 g 2 (d,p,x,e) 

with g±,g 2 T -periodic functions of x. Moreover, 

(i) If (8 e (x) , p e (x)) and (8(x), p(x)) are solutions of the original system (13. 2 p and averaged 
system (13. 5 p respectively, with the initial value such that 

| *(0)- #0)| + | 0.(0) -0(O)| = 0(e), 

then 

| p e {x) - p Q (x)\ + | 8 £ {x) - o (x)| = O(e), 

for times x of order 1/e. 

(ii) If Pq is an equilibrium point of (13.51) such that the corresponding Jacobian matrix 
has no eigenvalue equal to zero, then (13.21) admits a T-periodic solution (8 e (t), p £ (t)) such 
that \(p e (t),8 e (t)) — P \ = 0(e), for sufficiently small e; if Pq is a hyperbolic equilibrium 
point of (13. 5p . then there exists eo > such that, for all < e < Eq, system i3.S\) possesses 
a hyperbolic periodic orbits J e (x) = Pq + 0(e) of the same stability type as Pq. 
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(iii) If (9 e (x), p s (x)) G W s (j E (x)) is a solution of system ( EJ| ) lying in the stable manifold 
of the hyperbolic periodic orbit j e (x) = Pq + 0(e), (9(x), p(x)) G W s (Po) is a solution of 
system i3.5\) lying in the stable manifold of the hyperbolic equilibrium point Pq and 

\(9 £ (0)),Pe(0)-m,Pm\ = O(e), 

then 

\(9 £ (x),p £ (x))-(9(x),p(x))\=0(e), 

for x G [0, +00). Similar results apply to solutions lying in the instable manifold on the 
interval x G (— 00, 0]. 

(iv) If Pq is an equilibrium point of system (13. 5p and there exists a neighborhood U(r; Po) 
of Pq ( with radius r and center Pq ) such that there is not another equilibrium point in the 
closure of U and 

deg(F,f/,P )^0 

with F = (Ft, F 2 ). Then for \e\ > sufficiently small, there exist a T -periodic solution (p e (x) 
of system \3. || ) such that 



ip £ (-,e) — > P as e — > 0. 



Proof. The proof of (i)-(iii) follows directly from [21| or [22[. The proof of (iv) is based a 
framework of coincidence degree theory. Without loss of generality, we assume Pq = 0. We 
define homotpoy operator U : C([0,T],M 2 ) x [0, 1] ->■ Li([0,T],M 2 ) by 

H(9, p, A) := XeF(9, p) + (1 - X)e 2 G(9, p, x, e), 

where G := (gi,g 2 ). According to 0, Ch. VI] , H is L-compact on x [0,T], where f2 is a 
bounded open set of C([0, T], M 2 ) defined by 

fi:={(^p)GC([0,T],M 2 ):||(^p)||<r}. 

We remark that (9, p) G is a T-periodic solution of system (I3.6P if and only if (9, p) is 
a solution of C(9,p) = H(9,p,0) in Q. Since there is not another equilibrium point in the 
closure of U, we let 

M 1 = min \F 1 (9,p)-F 2 (9,p)\>0, 

(6,p)edu 

M 2 {e) = emax{|5fi(^,p,x,£:) - 5-1 (9, p, x, e) \ : (9,p,x,e) eU x[0,T}} 

with M 2 (e) — > as £ — )• 0. We also assume M 2 (e) > 0, otherwise Po is a solution of system 
(13. 6 p and the result is proved. 

First, we claim that, for each e G (— e ,0) U (0, e ) with e = max ee [_ ljl ]{M 1 /M 2 (e)}, 
there exists no solution (9, p) G dVt for the operator equation 

£(9,p)=H(9,p,\), AG (0,1]. (3.7) 



In fact, if (9,p) G dfl is a solution of ( 13. 7p . then there exists £ G [0, T] such that 

^'(0+p'(0 = o, 

|| (0, p) || = max J&(x) + p(x) = Je^)+p\i) 

x£[0,T] v v 

and 

6>(0p'(0 < o. 

Thus, it follows that 

o = |^(0 + ^(0l 

> |e| • \F X {9, p) - F 2 (9,p)\ -e 2 \ 9l (9,p,x,e) - g 2 (9,p,x,s)\ 

> \e\Mi - eM 2 {e) > 0, 

which is a contradiction. 

Without loss of generality, we suppose that 

£(9,p)^H(9,p,\), (9,p)edQ (3.8) 

holds for A G [0, 1]. Otherwise, the result is proved for (9,p) G dQ. Thus, we can apply the 
homotopy property of the coincidence degree and obtain 

\D C (C - H(; 0), Q)\ = \D C {£ - H(; 1), Q)| 
= |deg(F, fiHR 2 , 0)| ^ 0. 

Hence, by the existence property of the coincidence degree, there is (9, p) G fi such that 
C{9, p) = %(9, p, 0). Then (9, p) is a T-periodic solution of ( I3.6p . Thus, owing to the change 
of variables in this theorem, there is a T-periodic solution (9,p) for system (I3.5p . □ 

We remark that the proof of part (iv) of Theorem 13.11 does not need the smoothness 
condition upon F. So, it is convenient to deal with the existence of periodic solutions for 
nonlinear systems with loss of smoothness by Theorem 13.11 

4. Periodic orbits and stability 

To study the dynamics of MAWs or QMAWs for system (11.11) . we must investigate the 
behavior of the periodic orbits for system (13.11) including the existence and stability. Ac- 
cording to the method of averaging, the equilibrium point of the averaged system determines 
the properties of the periodic orbit of the corresponding perturbed system. For example, 
the equilibrium point with its eigenvalue of linearization nonzero implies that there exists at 
least one periodic orbit; in addition, if the equilibrium point is hyperbolic, then the periodic 
orbit has the same type of stability as the equilibrium point, for sufficiently small parameter 
e. 
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In order to find periodic orbits of system (I3.2p . it is sufficient to find equilibrium points 
of the averaged system (13. 5p . In the following, we will discuss the existence and stability of 
the equilibrium points for the averaged system (13. 5p . For simplification, we assume that 

2 =7T + 0i, 0i = -7r/2, #0 = ^4/3, v = -B/2. 

Recalling system (13.51) . together with the assumption, we have the averaged system 



< 



1 c\ m 2 ( -, (4.1) 



Notice that there exists a constant cq > such that for each c G (0, cq), the equation 

P 4 - 2§ p 2 + j = (4.2) 
has at least two real roots 



Moreover, equation (14.21) implies that 

A r 4 r 2 

^-|r)-W-^) = o. 

Thus, we can find four equilibrium points as follows 

n / , [FT [B~ 2 & , 7T 0x 



The eigenvalues of the equilibrium points Pi 2 and P3 ; 4 are given by 



' 2v^ V A 2 5 1,2 45p 2 2 

and Ag 1 ] = — A 12, A^ 2 ] = — A^, respectively. So, the equilibrium points Pi j2 and P 3j4 are 
hyperbolic, and as a consequence persist as periodic orbits for system 03.2p ; in addition, 
these periodic orbits are instable. 
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If 9 = &7r — (f>i/2, to find the equilibrium points, one will solve the following algebraic 
equation 



2 KI 5 2 ' 



W + A{? + -)- B?{? + -) = 0. 
Equation (14.31) has at least two roots 

P3 ,4 = ±^+o(c) e (-oo, -^V5) U (^75, +oo), 
for sufficiently small integration constant c > 0. 



■rid 



c 2 IS 



.* A 




4 2 



2" T 




(4.3) 



Figure 1: The phase portrait associated with the averaged system (|4.ip . All the equilibrium points persist 
as periodic orbits for system 



If 9 = kir + tt/2 — <pi/ 2, the algebraic equation 
needs to be solved. Note that 



/ ± 



4-§^ 3/2<o and ^ ±oo )= +oo > 

9 



(4.4) 



for sufficiently small positive constant c. By the mean value theorem, equation (14.41) has 
two roots p5 t 6 such that 



for sufficiently small c. As a consequence, four equilibrium points of the averaged system 
( 14. ip are obtained as follows 



The eigenvalues of the linearization at the equilibrium points P§£ and P^g are given by 



respectively. The equilibrium points P^s imply that that two instable periodic orbit of 
system ( 13. 2 p exist; while the equilibrium points P^q are nonlinear centers, and also persist 
as periodic orbits for system ( 13. 2p . The phase portrait for system ( 14. ip is given in Figure 
1. Since the periodic orbits corresponding to the equilibrium points P^q are not hyperbolic, 
one can not conclude their stability. This question is left open for further study. 

In summary, for the equilibrium Pj (i — 1, • • ■ , 4, 7, 8), there exists Eq > such that, for 
all < e < e , system (13. 2p possesses a unique hyperbolic periodic orbits %(x) = Pj + 0(e), 
which is instable. 

We also remark that, if (9 £ (x), p e (x)) G W s {^ £ (x)) is a solution of system ( 13 .2p lying in 
the stable manifold of the hyperbolic periodic orbit j E (x) = Pi + (D(e), (9(x), p(x)) e W s {Pq) 
is a solution of system ( 14. ip lying in the stable manifold of the hyperbolic equilibrium Pj (i = 



,4,7,8) and |(0 e (O),p e (O)) - (0(O),p)(O)| = 0(e), then \(9 E (x), p E (x)) - (6(x), p)(x)\ = 



0(e), for all x e [0,+oo). Similar results apply to solutions lying in the instable manifold 
on the interval x G (— oo, 0]. 

Although we study system (14. ip for sufficiently small c > 0, c also can be taken on a 
open interval (0,c), for some positive constant c. By continuous dependence of solutions 
with respect to the parameters, there is a connected set C of T-periodic solutions for system 
( 13. 2 p and then for system (13.11) . Since 

ip(t,x) = P(x)expz[6(x) — fit] 

= R(x) ( cos[6(x) + vx — pt] + i sin[6(x) + vx — pt}) , 
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and 




where 



1 




xo+T 



C 



and Q(x) = Q(x) — v is a T-periodic function with zero mean value, whether ip(t,x) is a 
MAW or QMAW depends on the choosing of the integration constant c. Precisely, If 2tt/i/ 
and T are rationally related, then ip(x, t) is a MAW; if 2tt/u and T are rationally irrelevant, 
then i/)(x, t) is not periodic but quasi-periodic, which is corresponding to a QMAW with the 
frequency u = (2tt/v,T). 

5. Numerical simulation 

To demonstrate the process of averaging to BECs, a specific example of numerical com- 
putation is given in the following. We take 

,,13 
9K X ) = 2 _ 2 Sm ' 

V(x) = 2sin2x - 1 

and the parameters 8 — l,c — 0.5, e = 0.01. Obviously, g(x) crosses zero. We can find 
equilibrium points for system (14. ip in the (6, p)-coordinates as follows 



Pjfi are nonlinear centers with eigenvalues of the linearization = Ag = He. 

Using the transformation these equilibrium points in the (R, 5')-coordinates with 
x = are given by 



P 5 , 6 : (±0.80, ±0.48), P 7j8 : (±0.98, =f0.82). 

We plot the solutions of system (13. ip starting from P ir i = 1,2, ••■ ,8, according to the 
averaged theorem, which are a good approximation to the periodic orbits, see Figure 2. 

6. Conclusion 

In conclusion, we have presented first-order averaging theorem in the periodic case for 
dynamics of MAWs (or QMAWs) in collisionally inhomogeneous BECs. The transformed 
system is non-Hamiltonian, and we indicate how the averaging theorems can be used to 
prove the existence and stability of periodic solutions. The questions as mentioned in the 
introduction have been answered. When the sufficiently small scattering length a(x) varies 



Pi }2 : (kir, ± 1.60), P 3A : (kn + -, ± 1.60) ; 

jIX 7T 

P 5 , 6 :(kn + —, ±1.02), P 7 , 8 : (Attt + -, ±1.33). 



P 1)2 : (±1.60,0), 



P 3 , 4 : (±0.31, 0) 
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Figure 2: The solutions portrait of system p.l| with initial value (R(Q), -R'(O)) = Pi, i = 1, 2, • • • , 8. Here, 
we take g(x) = \ — \ sin 2a;, V(x) = 2 sin 2a; — 1 and the parameters 8 = 1, c = 0.5, £ = 0.01. These solutions 
are a good approximation to the periodic orbits for x of order 1/e. The solutions with solid lines are instable. 



periodically in spatial variable x and crosses zero, infinitely many (positive measure set) 
MAWs and QMAWs can be proved to exist by adjusting the integration constant c on some 
open interval. 

A numerical approximation of periodic orbits is given for some prescribed parameters. 
We remark that, expanding at each equilibrium point and combining with multiple scale 
perturbed theory, such as work in 0, [24 26] , there may be a better approximation for each 
continuation periodic orbit. However, we emphasis on the theory frame of averaging to study 
dynamics of collisionally inhomogeneous BECs. 

In the end, it should be remember that the asymptotic approximation are valid for small 
e, but how small is usually a difficult problem. However, one advantage of averaging is 
obvious that it is set up for an easy return to the original variables. 
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